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Abstract 

We perform numerical simulations of vortex motion in a trapped Bose-Einstein condensate by 
solving the two-dimensional Gross-Pitaevskii Equation in the presence of a simple phenomenological 
model of interaction between the condensate and the finite temperature thermal cloud. At zero 
temperature, the trajectories of a single, off - centred vortex precessing in the condensate, and of a 
vortex - antivortex pair orbiting within the trap, excite acoustic emission. At finite temperatures 
the vortices move to the edge of the condensate and vanish. By fitting the finite -temperature 
trajectories, we relate the phenomenological damping parameter to the friction coefficients a and 
a , which are used to describe the interaction between quantised vortices and the normal fluid in 
superfluid helium. 

PACS numbers: 03.75.Lm, 47.37.+q 
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I. INTRODUCTION 



Atomic Bose - Einstein condensates provide an ideal testing ground to study quantised 
vorticity. The generation, dynamics and decay of simple vortex configurations have been 
observed and described l[ 2, [J • At finite temperatures, the interaction of the condensate 
with the thermal cloud damps the motion of structures (collective modes, solitons, vortices) 
5|, |6j. Current investigations of atomic Bose-Einstein condensates are concerned with this 



damping 
damping 



is y, q, 
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121 ]. Our aim is to use a very simple phenomenological model of 



131 ] and relate the value of the damping parameter 7 to the friction coefficients a 



and a', w" 



rich in the context of superfluid helium, are used to represent finite temperature 



energy 



effects [17|,L1EL19|,12Q|. It is well known, that even at absolute zero, whereas the tota' 
of the condensate is constant, the kinetic energy of a vortex structure can decrease |2l| due 
to the generation of sound. In assessing the finite - temperature damping, we shall therefore 
isolate the transformation of kinetic energy into sound energy. 



II. MODEL 



We solve numerically the following two-dimensional Gross-Pit aevskii Equation (GPE) 
which governes the time evolution of the order parameter, ip(r,t) = |^(r, t)|e* 5 ^ r '^: 



ifhfi = (-^ V2 +9\^\ 2 + Krap " V>, (1) 

where r is the position, t is the time, 

p(r,t) = |^(r,t)| 2 , (2) 

and 

v = -VS(v,t), (3) 
m 

are respectively the density and the velocity. The other quantities which appear in Eq. [1] 
are the coupling constant, g = 4irh 2 a/m, the atomic mass, m, the scattering length, a, 
the chemical potential, fi, and the trapping potential, V tr ap- To simplify our calculation, 
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hereafter we assume that the condensate is two dimensional and the trapping potential has 
the form 



Vtrap = 2 mUJ l + V 2 ) 



(4) 



where u± is the angular frequency of the trap In order to perform our analysis, we decompose 
the total energy, E tat , into kinetic, internal, quantum and trap contributions, 



where 



Etot — Ekin + Ei n t + Eg + Etrapi 



— (x/p(x,t)v(x,t)J d 2 r, 



(5) 



(6) 



Eint{t) = I 0(p(x,t))Vr, 



E q (t) = J ^(vv^M) 2 A, 

E t rap(t) = J p(-x,t)V tr d 2 r. 



(7) 



(8) 



(9) 



Furthermore, we decompose the kinetic energy, Euni into a part due to the sound field, 
E sound, and a part due to vortices, E vortex : 



E, 



kin 



E 



sound 



+ E, 



vortex • 



(10) 



At a given time t, the vortex energy, E vortex , is obtained by relaxing the GPE in imaginary 
time, which yields the lowest energy state for a given vortex configuration 22| . At this 
point the sound energy is recovered from E sound = E kin — E vortex . We rewrite the GPE in 
dimensionless form in terms of the harmonic oscillator energy hw±, the harmonic oscillator 
length y/hj (2rruJ±) and the harmonic oscillator time uj 1 . We obtain: 



\v 2 + v trap + c\ 



(11) 



where 
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v trv = \(^ + y 2 ) • (12) 

C-±2i, (13) 

N is the number of atoms and L the extension of the condensate in the z - direction. 
Unless stated otherwise, we set C = 2000 in our calculations. The added phenomenological 
dissipation parameter, 7, models the interaction of the condensate with the thermal cloud 



3, Ql 



Its microscopic justification was provided by Penckwitt et al 15] and Gardiner 



at al 16]; they studied the growth of the condensate in the presence of a rotating thermal 
cloud and found Eq. dH with 7 = 4mga 2 kT/(irh 2 ) = HW + /kT « 0.01, where W + is the 
rate at which thermal atoms enter the condensate due to collisions, k is the Boltzmann's 



constant and 5 = 3 a correction factor. Eq. [TTJ was also derived by Choi at al 11] and 



Tsubota et al [13j. The numerical method which we use is based on second order, centred 
finite differences in x and y, and Crank-Nicholson time stepping. The time stepping involves 
operator factorization of the right-hand side of Eq. [TTJ the Thomas algoritm, and an inner 
iteration to cope with the nonlinearity. The typical discretization is At = 10~ 3 in time and 
Ax = Ay = 0.087 in space; the convergence was checked by test runs with smaller values. 
The chemical potential, which can be initially estimated as 

from the Thomas-Fermi solution of the dimensionless GPE, is found by the single-step 
integration: 



In 



(\m\ 2 ) 

<|V(t+At)P) 



,m^r> (15) 

p 2At K ' 

after the initial Thomas-Fermi solution has relaxed to a time - independent solution in the 
harmonic trap, where (...) denotes spatial average. The calculation is performed in a square 
box of size D. We choose D so that it is larger than the trapped condensate, and impose 
boundary conditions ip = at x = ±D/2, y = ±D/2. Typically, D = 13. Fig. [1] shows a 
typical condensate without any vortices. 
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III. DECAY OF ONE VORTEX 



In order to imprint a vortex at location (xo,yo), we take for initial condition if) the 
Thomas- Fermi approximation, multiplied times a suitable function which is proportional to 
(x — x ) + i(y — dq) and vanishes at (x ,y ). Without dissipation (7 = 0), if the vortex is 
initially placed off-centre, it precesses around the trap with period 2it/ujq, following an orbit 
of constant total energy under the action of the Magnus force {t| , as shown in Fig. [21 The 
error in locating the vortex is of the order of ± 0.01. The angular frequency of the orbit u)$ 
increases with x as shown in Fig. [3j For example, if Xq = 0.9, y$ = 0, the precession is an 
orbit of total energy E tot ~ 17.42. The kinetic, quantum, internal and trap contributions 
to E tot oscillate with time; their average values are: E kin ~ 0.08, E q ~ 3.55 x 10~ 6 , E int ~ 
8.57 and E tra p — 8.77 respectively. The relative size of the oscillations during an orbit are 
respectively AE km /E kin ~ 0.125, AE q /E q ~ 0.097, AE int /E mt ~ 0.029 and AE trap /E trap ~ 
0.031, where AE kin , AE q , AE int and AE trap are the amplitudes of these oscillations. The 
likely reason of these oscillations is that, as the vortex precesses, it generates sound waves, 
which, unable to escape the trap, are reabsorbed by the vortex. This is the likely reason 
of the orbital wiggles apparent in Fig. [2] (the other reason is that the centre of mass of a 
condensate containing a relative large vortex hole oscillates). The total energy, Ef Q t is a 
constant of motion; the finite discretization of the numerical scheme is such that over the 
duration of the run, the relative change is thus a small value. AE tot /E tot ~ 0.012, rather 



In the presence of dissipation (7 7^ 0), the vortex loses energy, spirals outward toward 
the edge of the condensate and then vanishes, as shown in Fig. H] and Fig. [5J It is apparent 
that, increasing 7, the decay of the vortex becomes faster. 



In the vortex filament model of Schwarz [19j the motion of a quantised vortex in superfluid 
helium is determined by the balance of Magnus and drag forces. The resulting equation for 
the vortex position s = s(t) is Schwarz's equation: 



than zero. 




-T7 = v s df + as' x (v n - v self - v s ) - as x s' x (v„ - v self - v s ) . (16) 



where a and a are temperature-dependent friction coefficients, v n and v s are the externally 
applied normal fluid and superfluid velocities and s' is the unit tangent vector to the vortex 
at s. In two - dimensions s = (x,y) and s = z, where z is the unit vector along the z 
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axes. In our case, v se if is the precessing velocity of the vortex in the trap in the absence of 
thermal cloud. Assuming v s = and v n = (stationary thermal cloud) and using cylindrical 
coordinates (r,9), the solution to Schwarz's equation is: 

r(t) = r(0)e- aulo \ 0(t) = 6(0) + w„(l - a)t. (17) 

By fitting the calculated vortex position (at given value of 7) to Eqs. [171 we deduce the 
friction coefficients a and a . The results slightly depends on the initial position of the 
vortex because the condensate is not homogeneous near the edge. Fig. [6] shows that the 
deduced values of a is approximately constant for x < 2 (centre part of the condensate) 
and decreases more rapidly for xq > 2 (outer part of the condensate). For initial condition 
sufficiently close to the centre of the condensate, we find that a is proportional to the 
dissipation parameter 7, as shown in Fig. [7J The transverse friction coefficient, a , is much 
smaller than a, thus more difficult to determine. Fig. [8] shows that a' is approximately 
proportional to 7 only for small values of 7. 

IV. DECAY OF A VORTEX PAIR 

We repeat the calculation for a vortex - antivortex pair initially located at (=br , 0), first 
for 7 = and then at increasing values of 7. Initially, the pair moves across the condensate. 
When the pair approaches the edge of the condensate, the two vortices separate, and 
move back toward the opposite side of the condensate, thus making a closed orbit before 
returning to the initial position, as showed in Fig. [9l 

It is well known that, in the absence of dissipation, a vortex-antivortex pair set at (± xq, 
0) in an infinite homogeneous condensate moves with (dimensionless) translational speed 
^00 = l/(2xo) = 1/do, where do is the initial separation distance between the vortices. We 
define v pa i r the measured velocity of the vortex pair when it moves near the centre parallel 
to the y axis and compare v pair with in Fig. [10J As expected, the best agreement is for 
small do. 

We now turn our attention to the energy. For a vortex pair initially located at (± 
1.43,0), we find E tot ~ 17.47. Because of the numerical resolution, E tot is not constant, but 
varies of a typical amount ± 0.015 over t = 140, which corresponds to four orbits of the 



6 



pair. That indicates a relative accuracy of ± 0.09 % in conserving the total energy. We 
used different initial positions in the case of a vortex pair and in the case of a single vortex. 
The quantities E int , E trap , E kin and E q have average values E int ~ 8.6, E trap ~ 8.8, E kin 
~ 0.11 and E q ~ 1.93 xl0~ 5 with oscillation of maximum relative amplitude AE int / E int 
~ 0.035, AE trap /E trap ~ 0.045, AE km /E kin ~ 0.127 and AE q /E q ~ 0.078. Fig. El and 
Fig. [T2l show the time dependence of Ei nt and E trap . By propagating the GPi? in imaginary 
time, we also compute separately the contribution of E sound and E vo p^ to E kin , see Fig. [13J, 
Fig. [Hand Fig. [151 Fig. [151 shows that the correlation coefficient [25(, cc, between E sound 
and E vortex is negative, cc ~ -0.844 for initial separation c? = 2.86. Indeed, close inspection 
of the time dependence of E sound and E vortex reveals that when E sound is approaching its 
maximum, E vortex is approaching its minimum. 

In the presence of dissipation, the total energy, E tot , of the vortex - antivortex pair decreases 
with time, and so do the contributions E kin , E int , E q and E trap to E tot . The energy decay 
is faster if the vortices of the pair have initial smaller separation. For example, Fig. fT7l 
compares the decay of E k i n for two vortex pairs, one with initial separation do = 2.86 and 
one with do = 1.5. The latter is faster, as the pair spends more time near the edge of the 
condensate. 



V. DISCUSSION 

Our calculation shows that some features of the motion of vortices in a Bose-Einstein 
condensate can be modelled relatively well using Schwarz's vortex dynamics. In particular, 
we have been able to relate our phenomenological damping parameter, 7, to friction 
coefficients a and a . 

The natural question which arises is then the relation between 7 and the temperature 
ratio T/T c , where T c is the critical temperature. To answer this question we use results 
of unpublished preliminary investigations [?J using the Zaremba-Nikuni-Griffin (ZNG) 
finite-temperature theory [81], which show that, for C = 500 and for a single vortex initially 
located at (xo,yo) = (1.3,0), the effective friction coefficient is a ~ 0.0018 at T/T c = 0.15 
and a = 0.0025 at T/T c = 0.267. Similar results were found for (xo,yo) = (0.65,0). Setting 
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now C = 500, we rerun our calculations of single - vortex trajectories with the same initial 
condition (xq,2/o) = (1-3,0) and find that we need to set 7 = 0.044 to obtain the same value 



a 



0.0020 of ref 



7|, and 7 = 0.08 to obtain a = 0.0025. We conclude that 7 = 0.044 and 



0.08 correspond respectively to T/T c = 0.15 and 0.27. 



Final 
model 



Y, our small values of a are consistent with the ZNG theory [?| and with Berloff's 
|23| . Further work will attempt to derive Schwarz's Eq.[T6]from the Gross-Pitaevskii 



equation in terms of 7. 
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FIG. 1: (Color online): Equilibrium condensate with no vortices. The density contours corre- 
spond respectively to 17%, 33%, 50%, 60% and 83% of the maximum density at the centre of the 
condensate. 
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FIG. 2: (Color online): Trajectory of a single vortex initially located at xq = 1, yo = in the 
absence of dissipation (7 = 0). 
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FIG. 3: (Color online): Angular velocity, loq of a vortex which orbits the trap in the absence of 
dissipation as a function of initial position (xo,yo). 
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FIG. 4: (Color online): Trajectories of a single vortex initially located at (xo,yo) = (0.9,0) for 7 
= (dotted circle in the central region — 1 < x,y < +1); 7 = 0.01 (solid spiral); 7 = 0.07 (dashed 
spiral); 7 = 0.1 (outer dotted spiral). 
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FIG. 5: (Color online): Radius of trajectory versus time of a single vortex with initial position 
(xo,yo) = (0.9,0) for (reading from the horizontal line at the bottom of the figure to the steepest 
curved line) 7= ; 0.004 ; 0.008 ; 0.012 ; 0.016 ; 0.02 ; 0.024 ; 0.028 ; 0.032 ; 0.036 ; 0.04. 




FIG. 6: (Color online): Friction coefficient a (circles) and a (triangles) as a function of initial 
position (xo,yo) for 7 = 0.003. 
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FIG. 7: (Color online): Friction coefficient a for a single vortex with initial position (xo,yo) = 
(0.9,0) (triangles) and (xo, yo) = (2,0) (circles) as a function of 7. The linear fit for a is: a = C1+C27, 
where ci=0.007 and c 2 =5.092. 
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FIG. 8: (Color online): Friction coefficient a corresponding to Fig. 7. 




FIG. 9: (Color online): Trajectory of left vortex of vortex - antivortex pair initially located at 
(x ,yo) = (±1,0) for 7 = 0. 
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FIG. 10: (Color online): Vortex pair velocity as a function of initial separation do (circles), mea- 
sured at the centre of the trapped condensate for 7 = 0, compared to the vortex pair velocity in 
an infinite homogeneous condensate (dashed line). The classical velocity of a pair of point vortices 



was also studied by Jones and Roberts 24]]. For larger initial separation distance our results are 
similar to their results (for example for do = 3.5, v our = 0.2857 and v t y ie i r = 0.3). For smaller 
distance our results are different. (For do = 1.78, v our = 0.5618 and = 0.4). 
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FIG. 11: (Color online): Dimensionless internal energy, Ei n t{t) = J C (/o(x, t)) 2 d 2 r, for vortex - 
anti vortex pair for initial separation do = 2.86 and 7 = 0. 
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FIG. 12: (Color online): Dimensionless trap energy, E tra p{t) = J p(x, t)Vt rap d 2 r, corresponding to 
FIG. 11. 
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FIG. 13: (Color online): Dimensionless kinetic energy, Eki n {t) = f \ (\/p(x, i)v(x, i)J d 2 r, corre- 
sponding to FIG. 11. 
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FIG. 14: (Color online): Dimensionless sound energy corresponding to FIG. 11. 
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FIG. 15: (Color online): Dimensionless vortex energy corresponding to FIG. 11. 
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FIG. 16: (Color online): Correlation between dimensionless vortex energy and dimensionless sound 
energy for vortex - anti vortex pair for 

do = 2.86 and 7 = 0. The correlation coefficient is cc = —0.844. 
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FIG. 17: (Color online) Dimensionless kinetic energy vs time for vortex - anti vortex pairs for 7 
0.003. Left: with initial separation do = 2.86. Right: with initial separation do = 1.5. 
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